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Magnetic susceptibility and muon spin rotation (/^SR) experiments have been carried out to study 
the effect of structural disorder on the non-Fermi-liquid (NFL) behavior of the heavy-fermion al- 
loy Ce(Ruo.5Rho. 5)2812. Analysis of the bulk susceptibility in the framework of disorder-driven 
Griffiths-phase and Kondo-disorder models for NFL behavior yields relatively narrow distribu- 
tions of characteristic spin fluctuation energies, in agreement with ^SR linewidths that give the 
inhomogeneous spread in susceptibility. /iSR and NMR data both indicate that disorder explains 
the "nearly NFL" behavior observed above ~2 K, but does not dominate the NFL physics found 
at low temperatures and low magnetic fields. 

PACS numbers: 71.27.-Fa, 75.30.Mb, 76.60.Cq. 



I. INTRODUCTION 



A. Ce(Rui_3,Rha:)2Si2 



The discovery of non-Fermi-liquid (NFL) phenomena 
in strongly-correlated electron metals raises fundamen- 
tal questians about the elementary excitations of these 
systemsJdfl The interest in NFL behavior is in large part 
due to ptlie expected robustness of Landau's Fermi- liquid 
theory,ad according to which interactions between elec- 
trons that do not precipitate a phase transition should 
not change the Fermi-liquid nature of the low-lying exci- 
tationSpi Many /-electron heavy-fermion alloys are NFL 
metals.EJ Attempts to explain NFL behavior often invoke 
the notion of a quantum critical point (QCP) at zero 
temperature, the critical behavior of which extends to 
nonzero temperatures and modifies the thermodynamic 
and transport properties of the metal. 

Recent nuclear magnetic resonance (NMR) and muon 
spin rotation (/iSR) investigations of NFL alloys have 
yielded unambiguous evidence that in some of these jpa- 
terials disorder is a major factor in NFL behavior.afl In 
such cases NMR and /iSR spectra reflect broad distri- 
butions of the local magnetic susceptibility x(r,T), the 
high-susceptibility end of which arises from regions of the 
sample that do not exhibit Fermi-liquid paramagnetism. 
It is clearly of interest to determine susceptibility dis- 
tributions in a considerable number of NFL systems, in 
order to understand better the systematic interplay be- 
tween quantum criticality and disorder in NFL behavior. 



The Ce(Rui_2;Rb2;)2Si2 alloy system exhibits a num- 
ber of magnetic and nonmagnetic ground states as a 
result of strong electron correlation. Q The end com- 
pound CeRu2Si2 is a Fermi-liquid heavy-fermion com- 
pound with no evidence for magnetic ordering, whereas 
CeRh2Si2 undergoes an antiferromagnetic (AFM) tran- 
sition at a Neel temperature T/v = 35 K to a state 
of local-moment AFM order. With decreasing x T/v is 
suppressed, and vanishes for x ~ 0.55. A second re- 
gion of AFM order in the phase diagram appears for 
0.05 ^ X ^ 0.25; here the ordering is between itiner- 
ant rather than localized electrons. NFL behavicg has 
been established for x in the neighborhood of 0.5.0 Re- 
cent measurements of the electrical resistivity and nis 
netic susceptibility of Ce(Ruo.5Rlio. 5)2812 below 1 
have been interpreted in, terms of a quantum Griffiths- 
phase NFL mechanismlld for magnetic fields ^1 T and 
quantum spin glass behavior at higher fields. 

At temperatures ^ 2 K ^^Si NMR measurernents in 
an aligned powder sample of Ce(Ruo.5Rho.5)2Si2lij have 
shown that the local susceptibility is inhomogeneously 
distributed. The width Sx{T) of this susceptibility dis- 
tribution was founxLto be in good agreement with both, 
the Griffiths-phasellj and so-called "Kondo disorder" Q'Ej 
models of disorder-driven NFL behavior, which predict 
essentially the same 5x{T). In these theories a charac- 
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teristic energy A is inhomogeneously distributed in the 
sample. In the GrifRths-phase theory A is the tunnel- 
ing energy Et associated with a spin cluster, whereas 
in the Kondo-disorder model A is the Kondo temper- 
ature Tji of an individual spin. Each model yields a 
distribution function P(A) that can be used to calculate 
sample averages of experimental quantities. For example, 
the sample-average "bulk" susceptibility x{T) is given in 
terms of the local magnetic susceptibility x(T; A) by 



X{T)= I dAP(A)x(T;A). 



(1) 



For simplicity xC^l ^) is often taken to be of the Curie- 
Weiss form 



X(T;A) = 



3(r + A) 



(2) 



where PcS is the /-ion effective moment number. 

In Ref. |ll| the data were analyzed under the as- 
sumption that the disordered susceptibility inhomogene- 
ity is correlated only over short distances ["short range 
correlation" (SRC), or static susceptibility correlation 
length ^ lattice constant a]. The sensitivity of the 
NMR spectral width to comes about because a given 
^^Si nucleus is coupled to a limited number of neighbor- 
ing Ce spins. If the susceptibilities of these Ce neighbors 
are uncorrelated or only slightly correlated because is 
short, then the interactions are averaged somewhat. If on 
the other hand is much longer than the distance to the 
Ce neighbors ["long range correlation" (LRC), a], 
then the local coupling is not averaged and the (frac- 
tional) width of the frequency shift distributian-is the 
same as that of the susceptibility distributionliilO 

Agreement with disorder-driven models was obtained 
in the NMR experiments under the assumption of SRC. 
It was further assumed that the coupling could be char- 
acterized by an effective number of Ce ions coupled 
with a fixed interaction strength to a given ^^Si nucleus. 
For best agreement rioff ~ 5, which is reasonable crystal- 
lographically. It has been shown, however, that spectra 
from a second NMR mjjslB^iK or /iSR can be used to test 
the SRC assumption.y-By 

This article reports measurements of magnetic sus- 
ceptibility and /iSR spectra for temperatures greater 
than ^2 K in a high-quality single crystal of 
Ce(Ruo.5Rho.5)2Si2, .-which complement the previous 
NMR measurements O Consistent with the NMR study, 
we find that the mean and width of the susceptibil- 
ity distribution are in agreement with disorder-driven 
mechanisms. The predicted behavior is "nearly NFL" 
rather than NFL, i.e., the local susceptibility is Curie- 
Weiss-like with a distributed local Curie- Weiss temper- 
ature A(r)/A;B, but the distribution function P(A) van- 
ishes at A = 0. Then at low temperatures the system 
should revert to Fermi-liquid behavior in the absence of 
other mechanisms. 
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FIG. 1. Temperature dependence of c-axis bulk (sam- 
ple-average) magnetic susceptibility in single-crystal 
Ce(Ruo.5Rho. 5)2812. Curve: fits to GrifBths-phase and 
Kondo-disorder models (indistinguishable on this plot). 



The results of our susceptibility and transverse-field 
/iSR (TF-/iSR) experiments may be summarized as fol- 
lows: 

(1) The bulk susceptibility data can be fit to the 
Kondo-disorder and GrifRths-phase disorder-driven NFL 
models, thereby yielding the parameters of the distri- 
bution functions P(A) that characterize these models. 
These distributions are najrower than found in our pre- 
vious NMR investigationaB and possess little low-energy 
weight; the Griffiths-phase distribution, in particular, 
yields zero weight|-at A = rather than a Griffiths- 
McCoy singularity.ll3 

(2) The /iSR linewidths agree with these distribu- 
tions in the LRC limit, indicating a macroscopic dis- 
tance scale to the susceptibility inhomogeneity respon- 
sible for the spread in muon Larmor frequencies. The 
SRC-limit inhomogeneity found in the NMR measure- 
ments apparently reflects differences in sample prepara- 
tion. These differences preclude the independent test 
of the range of correlation (SRC vs LRC) discussed 
above. The amount of disorder derived from the disorder- 
driven models explains the nearly-NFL susceptibility of 
Ce(Ruo.5Rho.5)2Si2 above ~2 K, but is not capable of ac- 
counting fpij-the NFL behavior found at low temperatures 
and fieldsJa^y 



II. RESULTS 



A. Static Susceptibility 



The temperature dependence of the bulk magnetic sus- 
ceptibility Xc in a field of 0.1 T applied parallel to the c 
axis is given in Fig. |l|. These data were fit to the sample- 
average susceptibility x{T) given by Eq. (^ with distri- 
bution function P(A) from both the Griffiths-phase and 
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Parameter 




NMR, 






single crystal^ 


aligned powder^ 


Grifffiths 


A 


1.8 ±0.2 


0.88 


Phase 


coo (K) 


41 ± 5 


170 




Pctt 


3.01 ±0.06 




Kondo 


9 


0.15 ±0.003 


0.16 


Disorder 


Sg 


0.017 ±0.002 


0.021 




Ef (eV) 


1 (fixed) 








3.04 ±0.07 





TABLE I. Parameters obtained from fits of Grifiitlis-pliase 
and Kondo Disorder disorder-driven NFL models to bulk sus- 
ceptibility data. See text for definitions. 
''Data of Fig. 



^Data of Ref. 



Kondo-disorder-models. The Griffiths-phase picture 
{A = Et) yieldia 



P{Et) cx 



0, 



Et> loq, 



(3) 



where wq is a high-frequency cutoff, so that P{Et) di- 
verges as Et ^ for values of the nonuniversal ex- 
ponent A < 1 (Griffiths-McCoy singularities). In the 
Kondo-disorder model the Zener exchange coupling con- 
stant g — pj^ where p is the density of conduction- 
electron states at the Fermi energy and J is the 
conduction-electron-local-moment exchange interaction, 
is given a modest Gaussian distribution around a small 
average value.aO The local Kondo temperature 



K 



exp(-l/|g|) 



can then be widely distributed because of its exponential 
dependence on g. 

In both models the local susceptibility x(r; A) is taken 
to have the Curie- Weiss form of Eq. (||) . A fit of Eq. (|l|) to 
the bulk susceptibility then determines the parameters of 
the distribution function, which can be used to calculate 
the average of any function of A over the distribution. 
The fits to the Griffiths-phase and Kondo disorder models 
are given by the curve in Fig. |l|, and are indistinguishable 
on this plot. 

The fit parameters so obtained (A, cutoff energy 
and pcff for the Griffiths-phase fit; average g and standard 
deviation 8g of g, Ep^ and pcff for the Kondo-disorder 
fit) are given in Table H together with the corresponding 
values from fits to the susceptibility of the NMR sam- 
ple (Ref. |ll|). The distribution functions P(A) from the 
present results are given in Fig. ^. 

The most important feature of these results is that 
for both models P(A) has low weight for small A. For 
Kondo-disorder fits this result is rip qualitative agree- 
ment with the NMR-sample datajlJ which also yielded 
narrow distributions. But the best Griffiths-phase fit is 
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FIG. 2. Distribution functions -P(A) of characteric ener- 
gies A in single-crystal Ce(Ruo.5Rho. 5)2812, obtained from fits 
of Griffiths Phase (Ref. 0, solid curve) and "Kondo disorder" 
(Refs. I and |l|, dashed curve) disorder-driven NFL theories 
to bulk susceptibility data of Fig. |l[ 



obtained for A > 1, i.e., for a zero rather than a singu- 
larity at P(£'t=0) [cf. Eq. (||)], whereas the NMR-sample 
susceptibility was best fit with a singular form of P{Et) 
(A — 0.88). This indicates that within the framework 
of these disorder-driven models for NFL behavior the 
/iSR sample exhibits only "nearly-NFL" behavior, since 
nonzero values of P(A=0) are necessary if the lowest- 
lying excitations are to be local-moment-like or cluster- 
like rather than Fermi-liquid in character. The theoreti- 
cal distributions are simply too narrow io. yield the NFL 
behavior observed at low temperatures .Ba 



B. NMR and /xSR Frequency Shifts 



In a paramagnetic sample spin probe (nucleus or 
muon) frequency shifts reflect the electronic spin polar- 
ization to which the spin probes are coupled via a com-, 
bination of dipolar and indirect hyperfine interactions.1151 
A given spin probe (nucleus or muon) i is coupled to 
neighboring /-electron moments j, resulting in a linear 
relation between the spin-probe-fceouency shift Ki and 
the /-electron susceptibility Xj :liiHi3'EJ 



E 



where the coupling constants aij can be expressed in 
terms of corresponding coupling fields i/™"'' — N pBOHj] 
N is Avogadro's number if the Xj are expressed in molar 
units. We can write Hf°^^ as the sum of transferred- 

hyperfine and dipolar contributions H^j^^ and H^^^^, re- 
spectively. 

As described in Sect. | above and elsewhere,E: 
the relation between an inhomogeneous distribution of 
XjS and the resulting distribution of KiS depends on the 
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range of correlation of the inhomogeneous susceptibil- 
ity. In the two extreme limits of SRC |aad,LIlC the rms 
spreads SK and Sx are simply related:lliHl3llJ 

Sx - SK/a* , 

where the effective coupling constant a* depends on the 
correlation length 




(LRC) , 
(SRC) . 

In the LRC limit this yields 

{Sx/xkRC = SK/\K\ , 
and in the SRC limit we can write 

^x/XsRC = {aLRc/asRc){SK/\K\) . 



(4) 



(5) 



(6) 



We will use these relations in the analysis of our /iSR 
results. 
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FIG. 3. Dependence of sample-average muon c-axis 
frequency shift Kc on c-axis bulk susceptibility Xc i^i 
Ce(Ruo.5Rho. 5)2812, with temperature an implicit parameter. 
Circles: raw data. Triangles: data corrected for Lorentz and 
demagnetizing fields. The error bars include uncertainty in 
the value of the demagnetization coefficient. 



C. Transverse-Field /iSR Spectra 

The /iSR experiments were carried out at the M20 
muon channel at TRIUMF, Vancouver, Canada. A mag- 
netic field Hq was applied parallel to the c axis of the 
single-crystal sample, and the muon spin was oriented 
perpendicular to Hq. In this TF-/iSR configuration the 
width of the /iSR resonance reflects the distribution of 
muon frequency shifts in the sample as long as lifetime 
broadening (spin-lattice relaxation) is negligible, which 
has been confirmed by zero-field fiSH measurements. 113 

The spectra were fit to both Lorentzian and Gaussian 
lineshapes; values of the goodness-of-fit parameter x^ 
were not appreciably different for these choices. Figure ^ 
gives the dependence of the average (line centroid) frac- 
tional frequency shift Kc on the bulk c-axis susceptibil- 
ity Xc iri Ce(Ruo.5Rho.5)2Si2, with temperature an im- 
plicit parameter, before and after correction for Lorentz 
and demagnetizing fields. The curvature of the corrected 
Kc vs Xc could be due to a number of causes, including 
crystal-field effects and a Curie "taiill, in the susceptibil- 
ity from trace metallurgical phases.Eil Above ~20 K Kc 
is approximately proportional to Xc, with a slope that 
yields a coupling field 



(7) 



-0.71 ±0.06 kOe/nB- 



This is considerably more negative than j-the value 
—0.38 kOe/fj,B found in undoped CeRu2Si2ii3 and in- 
dicates that H^^^ is weaker in Ce(Ruo.5Rho.5)2Si2. 
This may refiect a difference in Ce-Si hybridization 
or, alternatively, a different muon stopping site in 



Ce(Ruo.5Rho.5)2Si2.t3 Indeed, different muon stopping 
sites are found in the end compounds CeRu2Si2 ((^ ^ 0), 
Wyckoff notation 26)E3 and CeRh2Si2 (two sites: (5 0), 
4c, and i^jO), 4d)@ 

The ratio of the fractional muon linewidth SKc to the 
frequency shift magnitude \Kc\ is plotted vs Xc in Fig. |^. 
The increase of SKc/\Kc\ with increasing susceptibility 
is a salieat-i'eature of disorder-driven models for NFL 
behavior Jj'tJ The curves in Fig. ^ give the fractional rms 
width Sx/x of the susceptibility distribution from the 



disorder-driven NFL distributions discussed in Sect. [I A 
where dx is obtained from the relation 



Sx{T)^\ / dAP{A)[x{T;A)-x{T)f 



1/2 



[Cf. Eq. i)]. ^ 

The data are in reasonable agreement with the 
disorder-driven pictures in the LRC limit, for which 
SK/\K\ (data points in Fig. ^ is the estimator of 
6x/x[Eq- (^)]. In the SRC limit this estimator is ob- 
tained by multiplying the experimental values of SK/\K\ 
by_aLRc/asRC [Eq. (|)]. The magniiudfij-of the slope 
of Kc vs Xc gives the value of aLRc ■tilullj It is more 
difficult to obtain agRc, but we can estimate separately 
the dipolar and transferred hyperfine contributions to 
this quantity and hence to flLRc/asRC- For dipolar cou- 
pling the value of this factor obtained from lattice sums 
[Eq. (^)] is 1.51. For a constant transferred hyperfine 
coupling strength to rieff near neighbors the factor is 
y^Ticff — 2-2.2 foe. the crystallographically reasonable 
range Ucs — 4-5 .EJ Thus in the SRC limit the experi- 
mental estimator of Sx/x is significantly increased, and 
agreement with the theoretical curves is worsened, for 
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FIG. 4. Data points: dependence of estimator 8Kc/\Kc\ of 
fractional c-axis susceptibility spread &Xc/Xc c-axis bulk 
susceptibility x^, with temperature an implicit parameter, 
in Ce(Ruo.5Rlio. 5)2812. Curves: dependence of 5x/x on x 
from fits of Griffiths-phase (solid curve) and Kondo-disorder 
(dashed curve) models to bulk susceptibility. 

both dipolar and transferred hyperfine coupling. This 
is evidence that in the present single-crystal sample the 
LRC limit is appropriate, i.e., that there is little effect of 
atomic-scale disorder on the local static susceptibility. 



likely to be related to the NFL behavior observed in 
Ce(Rui_:rRha;)2Si2 at tempftiatures below 1 K and mag- 
netic fields below 1 kOe.ffl The latter exists over a 
range of Rh concentration x =p, 0.5-0.6, as expected 
from the Griffiths-phase picture,El but the large clusters 
with low characteristic energies needed in this scenario 
are simply not observed in magnetic resonance experi- 
ments at higher temperatures and fields. Furthermore, 
it is unlikely that Grifhths-phase clusters form only for 
T < 1 K, since this temperature is considerably smaller 
than the energy scale for the disorder (~30 K, cf. Fig. ||). 
Magnetic-resonance shift and linewidth measurements 
would be desirable as direct probes of disorder in the 
low-temperature low-field NFL region, but unfortunately 
accurate data cannot be obtained at the low fields ('-^lOO 
Oe) required for the NFL behavior. 
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III. CONCLUSIONS 



The fractional width 5xc/Xc of the inhomogeneous sus- 
ceptibility distribution in Ce(Ruo.5Rho.5)2Si2 is consid- 
erable (~100% at low temperatures, see Fig. ^), but the 
parameters of the distribution functions (Table |) indi- 
cate smaller widths (smaller c^o and larger A > 1 for the 
GrifRths-phase fit; smaller 5g/'g for the Kondo-disorder 
fit) in the /iSR single crystal than in the NMR aligned 
powder sample of Ref. The origin of the disorder 
may also be different in the two samples, since the fits 
agree better with the LRC limit in the present work but 
with the SRC limit in the NMR study. This suggests that 
local disorder in the Ce-moment-conduction-electron hy- 
bridization strength is significant in the NMR sample but 
relatively weak in the /iSR single crystal, although varia- 
tions over length scales longer than a few lattice param- 
eters remain as a source of inhomogeneity in the latter 
sample. 

In both samples the low probability of small en- 
ergy scales (cf. Fig. ^ is suggestive of low spectral 
densities of spin fluctuations, quantum or thermal, at 
low frequencies. This is in agreement with the rela- 
tively slow /iSR spiuJattice relaxation rates observed in 
Ce(Ruo.5Rho.5)2Si2,tII since the muon relaxation rate is 
proportional to the strength of the thermal noise spec- 
trum at the low muon Larmor frequency.E3 

The data taken as a whole indicate that the disor- 
dered susceptibility revealed by the magnetic resonance 
(/iSR and NMR) experiments, while significant, is not 
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